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1 Introduction 

The Gross-Pitaevskii equation (GPE), derived independently by Gross [9^ and Pitaevskii 
K*" I [18], arises in various models of nonlinear physical phenomena. This is a Schrodinger- 

■ type equation with an external field potential Vext{t, x) and a local cubic nonlinearity: 

(N ■ ^2 ^ 

'sj" ' ihdfU^ l^U — VextU + (3 \u\ u. (1.1) 

• ' 2m 



o 



Abstract 

In this paper, we establish the global well-posedness of the Cauchy problem 
for the Gross-Pitaevskii equation with an angular momentum rotational term in 
which the angular velocity is equal to the isotropic trapping frequency in the space 



The GPE (|l.ip in physical dimensions (2 and 3 dimensions) is used in the mean- 
field quantum theory of Bosc-Einstcin condensate (BEG) formed by ultracold bosonic 
coherent atomic ensembles. A rigorous derivation of the GPE for the dynamics of (non- 
rotating) BEG has been obtained by Erdos, Schlcin and Yau 7, 8J. Recently, several 
k>( \ research groups [TOl [151 [HI [17] have produced quantized vortices in trapped BEGs, and 

5_j ■ a typical method they used is to impose a laser beam on the magnetic trap to create a 

I harmonic anisotropic rotating trapping potential. Seiringer has discussed the stationary 

GPE in [20 , Lieb and Seiringer have rigorously derived for the description of the ground 
state asymptotics of rotating Bose gases in (TH] . The properties of BEG in a rotational 
frame at temperature T being much smaller than the critical condensation temperature 
Tc |12j are well described by the macroscopic wave function u{t,x), whose evolution is 
governed by a self-consistent, mean field nonlinear Schrodinger equation (NLS) in a ro- 
tational frame, also known as the Gross-Pitaevskii equation with an angular momentum 
rotation term: 

ihdtu H Aw = V(x)u + NUq IuP u - ni^u, x € R^, i ^ 0, (1.2) 

2 m 

where the wave function u{t, x) corresponds to a condensate state, m is the atomic mass, 
h is the Planck constant, N is the number of atoms in the condensate, is the angular 
velocity of the rotating laser beam, and V{x) is an external trapping potential. When 
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a harmonic trap potential is concerned, V{x) — ^ ^j^'j with uji, uj2 and being 

the trap frequencies in the xi-, X2- and xs-direction, respectively. The local nonlinearity 
term NUo \u\^ u arises from an assumption about the delta-shape interatomic potential, 
where Uq — Airh'^as/'m describes the interaction between atoms in the condensate with 
Os (positive for repulsive interaction and negative for attractive interaction) the s-wave 
scattering length. Lz = ~'ih{xidx2 ~ X2dxi) is the third component of the angular 
momentum L ^ x x P with the momentum operator P — ~ih\7 . 

After normalization, proper nondimensionalization and dimension reduction in cer- 
tain limiting trapping frequency regime [19 , the system (|1.2p becomes to the dimen- 
sionless GPE in d-dimensions {d = 2, 3): 

1 9 ^ 

iut + -Au ^ Vd{x)u + Pd \u\ u - nizu, X eM.'^,t>0, (1.3) 
where — i{xidx2 — X2dxi) and 



with 7i > 0, 72 > and 73 > constants, /3 = , uq = y and uj,n = 

min{wi, W2, W3}. 

In general, it is a rather complicated process about the dynamics of solutions (in 
particular, vortex) for GPE (|1.2p under the interaction of trapping frequencies and 
angular rotating motion. The recent numerical simulation of GPE ()1.2p for different 
choice of trap frequencies (71, 72) can help us to understand the complicated dynamical 
phenomena caused by the angular rotating and spatial high frequency motion. The 
case of different frequency 71 7^ 72 gives much complicated behavior and thus is rather 
difficult to be studied rigorously To our knowledge, the equation (|1.2p has been 

only investigated for some specific cases by numerical simulation. Therefore, to develop 
methods for constructing analytical solutions to the GPE p.ip or some specific cases is 
the first step in order to understand the dynamics caused by the trapping and rotation. 

To begin with, we first consider the case 71 = 72 = 73 = w which means the spatial 
isotropic motion. In order to derive the exact analytic formula for the solution to the 
linear equation, we have to assume that the angular velocity is equal to the isotropic 
trapping frequency, i.e., fl = to. In the present paper, we focus on the Cauchy problem 
of the Gross-Pitaevskii equation with an angular momentum rotational term in three 
dimensions 

1 CJ^ 

iut + -Au = — |x|^ u + (3 \uf u - ujLzU, x eR^,t^O, (1.5) 
u(0,x) = uo(x), X e (1.6) 

where the wave function u = u{t, x) : [0, 00) x M'^ — > C corresponds to a condensate 
state, A is the Laplace operator on R'^, w ^ 1 and [3 > are constants, and = 
—i{xidx2 —X2dxi) = i{x2dxi —Xidx2) is the dimensionless angular momentum rotational 
term. 

We assume that the initial value 

uo(x) e S :== {u e H^{R^) : \x\u£ ^^(R^)}, (1.7) 

with the norm 



\\u\\^^\\u\\^. + \\\x\u\\^2. 

Note that for multi-dimensional GPE ()1.2p . nothing is known about the exact in- 
tegration except for the case 71 = 72 {— 73 for 3D) considered in [H [S] without the 
angular momentum rotational term, namely, 17 = 0. 

In the case of two dimensions [TT], the linear operator idt + ^(A — |x|^ + 2ujLz) 
can be written as the form idt -t- ^(V — 16)^ where h = {—UJX2, luxx) satisfies the Coulomb 
gauge condition V-6 = 0. However, in three dimensions, we can not find such a potential 
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b satisfying the Coulomb gauge condition V • 6 = 0. So that it is impossible to write the 
linear operator as a similar operator as in the case of two dimensions, and this makes 
the problem be much more difficult than the two dimensional case. 

In addition, there are three ingredients that play important roles in the proof of our 
result. The first involves the solution of the Cauchy problem to the linear equation 

1 uj"^ 

iut + -Au = —\x\ u- LuL^u, X eR^, t ^ 0, 

u{0,x) = uo{x), X e R^, 

which is significant for investigating the properties of the evolution operator correspond- 
ing to the hnear operator idt + — ^ \x\'^ + ujLz- The second one is to obtain the 
Strichartz estimates for the foregoing linear operator. The last one is that there exist 
two Galilean operators J{t) and H{t) (as blow) which can commute approximatively 
with the linear operator and can be viewed as the substitute of V and x respectively in 
the non-potential case. 

Now we state the main result of this paper. 

Theorem 1.1 Let uq G S and p £ (2,6). Then, there exists a unique solution u(t,x) 
to the Cauchy problem p.5p - (|1.6p . And the solution satisfies, for any T £ (0,oo), that 

u{t, x), J{t)u{t, x),H{t)u{t, x) e C(M; L^{m.^)) L'^^"' (0, T; LP{M.^)), 

where ^ i ^ f ' '-^i^) '^'^'^ ^(^) '^''^ defined as below as in (j2.14p and (j2.15p . respec- 
tively. 

Remark 1.2 Since the GPE (jl.Sp (or p.3p ) in a rotational frame is time reversible 
and time transverse invariant, the above result is also valid for the case when t < 0. 

The paper is organized as follows. In Sec. O the evolution operator of the linear 
equation and the Strichartz estimates about the former operator are first established. 
Sec. |3] is devoted to the derivation of some conservation identities such as the mass, 
the energy, the angular momentum expectation, and the pseudo-conformal conservation 
laws in the whole space M.^ for (|1.5p - (|1.6p . Finally, the nonlinear estimates and the proof 
of Theorem II . II are obtained in SecjH 



2 The Strichartz estimates and some main operators 



On the analogy of the Schrodinger operators with magnetic fields in [T] [6] , we can also 
define the propagator associated with the self-adjoint operator 



-A + oj^ \x\'^ - 2ujL, 



i.e., S{t) = e *^ot^ which can be explicitly expressed as 

/ io \^ f iiol ^^^J^^ cot(Lut) — x-y 



Hy)dy (2.1) 

giw^cot(a)t) I ^iujA{t)x-v^iuj^-cot{ujt)^^y^^y^ {2.2) 

Uiy)dy, (2.3) 



27ri sin(wt) 

UJ 

2Tri sui{ujt) 

UJ 

2'Ki snY{Ljt) 

for < ^ ^ ^ through a complicated computation, where x := (— 2:2, xi, (csc(a;t) 
coi{ijjt))x3) , the matrix A{t) is defined by 

cot{ujt) -1 
A[t) = I 1 cot(cjt) 

csc(wi) 
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and A'^{t) is the transpose of the matrix A{t). Note that this formula is vahd only for 
small time, due to the singularity formation for the fundamental solution. 
Next, we derive the dual operator S*{t) of S{t). From 



{S{t)f,g) ^ {S{t)f{x),-g{x)) = {S{t)f{0,9i-0) 



sin(wt) 



|(Ar(t))-l({-„)| 



sin(a;t) 



sin(wi) 



iv) 



|(AT(t))-lc| 

e'-^ 2;^ '=°*('^*)/((A^(t))"iC/w) 



i sin(wi) 

U.S*{t)g) 



-.zcj-^^— cot{ujt) ^ { ^iuj^^ cot{ujt) 



g(-) md^ 



we can define the dual operator as 

ioj cot(cji) / —iujA'^{t)x-y iuj cot(cjt) 



S*{t)g:^ 



2'Ki sin(a;i) 



UJ 



2ni sin(ajt) 



g{-y)dy 



giw-l^ cot(wt) / ^iuiA'^(t)x-y^iui^^cot(uit)g^y-^^y^ ^2.4) 



In order to obtain the Strichartz estimates, we have to estimate the norm \\S{t)S* {s)\\ j^i^j^^ 
Indeed, we have the following proposition: 

Proposition 2.1 Let u{t) = S{t)uo, then u{t) satisfies the linear equation p.Sp . And 
the operator S{t) has following properties: 

(1) Sit) is unitary on L^ , i.e. ||'5'(i)||^2_,^2 = 1; 

(2) \\S{t)S*{s)\\^,^^^ < |i - s\-^'^ forQ<s<t^^. 



Proof. By computation, it is easy to see that 



lUt 



2'Ki sai{u}t) J Jjjs 



- cot{Lot) — x-y 



-iuj cot(ajt) 



+ u)'^ csc^ {ujt) J^^ — I- w^(csc(wi) — cot(a;t)) csc(a;i)a;3?/3) uo{y)dy, (2.5) 



1a / ^ 
-Am = — - 

2 \27ri sin([ji) 



-iuj cot{ujt) 



cot^(wi)^^-— ^ - + cot(u;t) x ■ y uo{y)dy, (2.6) 
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^ I |2 

\x\ u 

2 ' ' 



2Tii sin(a;t) 



— \x\ ) uo{y)dy, (2.7) 



and 



^ 27ri sin(wt) 
Summing (|2?5|l -([2?8 ll . we have 



- cot{ujt) — x-y 



{x2yi - 112/2) cot{ujt) 



+ xiyi + X2J/2 



uo{y)dy. (2.8) 



1 

TO* + 2~ ^ ^^^^ ~ 

which yields the desired result. 

For the unitariness of S{t) on i^, it is trivial by using the Plancherel theorem. Thus, 
we omit the detailed proof. 

Next, we prove the dispersive property of S{t). Indeed, we have for < s < i < ^ 



S{t)S*{s)f 



'(27r)^ sin(a;s) sin(cL;t) 



iu 2 cot{(jji) I iujA{t)x-y iuj ^ cot{ujt) 



E3 



^2 cot(a;s) / ^^t^-^ (s)y-2^ iLJ ^ ^ cot(ct. 



'^f{z)dzdy 



e 2 



-27r sin(a;s) sin(u;t) ^ 

•^-1 (^e^-^^cot(.a)^(.)^ (u:A^{s)y)dy 



e 2 



sin(a;s) 



-27rsin(ciji) 



|(^^(^.))-i5|' 



[cot(L^t)+cot(L^s)] ^-1 I gjw^cot(a)s) yj- -J j (^)d^ 

|(Ar(a))-l5|^ 



- [cot(Ljf) + COt(LJs)] 



— 27rsin(ct;t)^ 



— sin(ti;i) 
sin(a;s) 



- [cot [ujt)+ cot {ujs)] 



■-°^i^^)f{-:,))\(B{t,s)x) 



— sin(ti;i) 

Ipi-i ^ i— [cot([^t)+cot(ii)s)] 



(i?(<,s)x-y)e^-^-*(-^)/(-2/)dy 



where the matrix B{t, s) is given by 
B{t,x) = 



Bii{t,s) Bi2{t,s) 
-Bi2(i,s) -Bii(t,s) 

csc(ijj/:) sin(a;s) 



with 



Bii{t, s) — sin^(tJs)(cot(wt) cot(cjs) + 1), Bi2{t, s) — sin'^(a;s)(cot(aji) — cot(tjs)). 
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Noticing that 

|(AT(3))-ieP 

r-lJ- 2Z L[cot((^t)+COt(, 



--i2n)- 
- csc^ [ujs) 



^ — [cot(cjt)4-cot(a;s)] 



< 



Tr[cot{u}t) + cot(ws)] 

sin^ {ujs)[cot{ujt) + cot(cjs)] 
we can get, for < s < i ^ that 



, ^COt(ljt)+COt(lJs) 



< 



TT sin(a;<) sin(cL's) [cot(a;t) + cot(tt's)] / Jj^s 



ii/iIl. <ii+^r'/'ii/iiLi 



\f{y)\dy 



7rsin(w(t + s)) 

<ii-.^r'/'ii/iiL. 



since |sint| ^ f |<| for \t\ s: f . □ 
Thus, we can obtain similar Strichartz estimates to the hnear Schrodinger operator 
by the standard methods (c.f. iJJj) provided that only finite time intervals are 
involved (c.f. [50. 

Proposition 2.2 Let I be an interval contained in [0,Tr/Auj]. Then, it holds that 

(1) For any admissible pair {'^{p),p) (that is, 2/j{p) = 3(1/2 — 1/p) for 2 ^ p < 6), 
there exists Cp such that for any (f) (E L'^ 

\\Sm\\L-ri.Hl;L^)^Cpm^.. (2.9) 

(2) For any admissible pairs (7(^1), Pi) and (7(^2), P2); there exists Cp-^^p^ such that 



S{t)S*{s)F{s)ds 



Lt(pi)(/:Lpi) 



• in{s<t} 

The above constants are independent of I d [0,7r/4w]. 
The integral equation reads 

u{t) = S{t)uo -if3 f S{t)S* (s) \u\^ u{s)ds 



Pl,P2 l|-^llL-f(P2)' (/;LP2) ■ (2-10) 



(2.11) 
Jo 

Since the initial data belong to S, we naturally need the estimates of \/S{t)(p and 
xS{t)4>. In fact, from ()2.ip . we can compute and obtain that 

WS{t)(j) — iijjxcot{ujt)S{t)<t) — iujS{t){x cot{u!t) — x)(f>, 

where x — (—xi, x\, (cotwt — c&c ut^x-i) , and 

S{t)Vcl} = iu;{xcotiujt) + x)S{t)(j) - iio coi{ijjt)S{t){x4>), 

which yield 



\JS{t)(j) =cos{ujt)S{t)(cos(u}t)\J - sin(tji)V)0 

— iijjs\n.{ujt)S{t) [(cos(a;t)a: — sin(a;t)a;) 
xS{t)(j) =cos{ijjt)S{t) [(cos(a;t)a; — s\n{ujt)x)(j)\ 



(2.12) 
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— sm{ojt)S{t){cos{ujt)V — sm{ujt)V)4>, 

OJ 
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(2.13) 



where V = {—dx2,dxi, (cotwf — esc 101)8x3)- 
Thus, we have 



S{t){-iV)(j) 

u> sm{u!t){cos{(jt)x + sin{(jt)x) — i cos{u)t){cos{ujt)V + sin(a;t)V) S{t)(j), 



and 



S{t)ui!X(j) 

u) cos{ut){cos{u)t)x + sm{u)t)x) + i sm{u)t){cos{u)t)V + sm{u)t)'V) S{t)(j). 

For convenience, we denote 

J{t) = uj sm{Ljt){coa{u}t)x + sin(wt)x) — i cos(a;<)(cos(wt)V + sin(a;t)V), (2.14) 
and the corresponding "orthogonal" operator 

H{t) = u> cos{u>t){cos{(jt)x + sin(a;t);r) + i sin(a;f)(cos(a;t)V + sin(a;t)V), (2-15) 

which will appear in the pseudo-conformal conservation law and play a crucial role in 

the nonlinear estimates. 

By computation, we can obtain the following commutation relation 



1 

J{t),idt + ^ \xf +ioL^ 

H{t),idt + ^A-^\x\^+u;L, 



= Oj{t), 
= Onit), 



(2.16) 



where 



and 



Oj{t) = {0,0,2iLOsm{u!t){usm{wt)x3 — icos{ut)dx3)) , 



Onit) = {0,0,2iuism{wt){u)cos{ujt)x3 + i sm{uit)dx3)) 



It is clear that for < t ^ 



\Oj{t)u\ 



^ „ 2iu! sm(u;t) ^ , , 
0,0, ^J z{t)u 



2 cos(a;i) — 1 



2lo sin(a;t) 



2 cos(wi) — 1 



Jz(t)u 



2u?t 
V2-I 



\J3{t)u\<UjH\J{t)u\ 



and 



\OH{t)u\ = 



0,0,- — ; \ ii^{t)u 

2 cosyujt) — 1 



<u?t\H{t)u\ 



(2.17) 



(2.18) 



where J2,{t) and H-iif) are the third component of the operators J{t) and -ff(f), respec- 
tively. 



In addition, denote M{t) = g-^'^T-t^^M) and Q{t) = e''^^'=°^('^*\ then 

J{t) = -icos{ujt)M{t){cos{Lot)V + sm{ujt)V)M{-t), 
H{t) =ism{ut)Q{t){cos{ujt)V + sm{(jjt)V)Q{-t). 



(2.19) 



Hao, Hsiao, and Li 



3 The conserved quantities 



Proposition 3.1 Let u be a solution of the equation (|1.5p with the initial data <j) G 
E(R'^). Then, we have the following conserved quantities for all t ^ 0; 
(1) The L^-norm: 



\\um, = \\uo\\,. 

(2) The energy for the non-rotating part: 

1 2 Lu'^ 2 4 

Eoiu) = - ||Vu||2 + Y ||a;u||2 + - \\u\\^ = E^oK)- 

(3) The angular momentum expectation: 

{L,){t) = f uL.udx ^ {L,){0). 

(4) The pseudo-conformal conservation law: 

\\Jit)u\\l + \\H{t)u\\l + Acos{ujt){l - cos{ujt)) II 
+ P\\u\\l = 2Eo{uo). 
Proof. For convenience, we introduce 

1 uj'^ 2 2 

eq{u) :— iut + -An ^ l-^l u — f]\u\ u + ujLzU. 

It is clear tlrat p.ip holds by applying the L^-inner product between eq{u) and u, 
and then taking the imaginary part of the resulting equation. 

Since wc can use the identity (j3.3p in the proof of p.2p . we derive (|3.3p first. Differ- 
entiating {Lz){t) with respect to t, and integrating by parts, we have 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



dt 



[ut{x2dxiU - xidx^u) + u{x2dx^ut - xidx^ut)] Ax 
-iut{x2dxiU - xidx2u) - iut{x2da:j^u - xidx^u)] dx 



1 OJ 

-Au — — \xf u- 13 \u\^ u + wL^u ) {x2dxiU - xid^^u) 
1 

-Aw — \x\^ u- 13 \u\^ u + ojLzU ) {x2dxiU - xidx^u) 



dx 



Re{Au{x2dxiU - xidx^u)) - -y (l^^l^ {x2dxt \u\^ - xid^^ |w|^)) 
- /3Re(|u|^ {x2dxi \u\^ - xid^^ |u|^))da; 



/ {2x1X2 |w|^ - 2x2Xi |u|^))da; 



=0, 

which yields the desired identity 

Next, we prove the energy conservation for the non-rotating part p.2p . We consider 

Re(e(7(w), Ut) = 0, 

where (•, •) denotes the L^-inner product. From the above, we can get 



/ 



|V-u|^ + '^dt |xu|^ + ^dt \uf + '^dtiuL.u) 



dx ^ 0, 



which implies the identity p.2p with the help of p.3p . 

Finally, the pseudo-conformal conservation law p.4p can be easily derived from the 
definition ()2.14p and (|2.15p of the operators J{t) and L[{t) with the help of the energy 
conservation for the non-rotating part p.2p . We omit the details. □ 
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4 Nonlinear estimates and the proof of Theorem 11.11 

With the help of dHHl), we can get 



j{t) |ur It = 2 |ur j{t)u - u^j{t)u, 

which imphes, in view of^+£=| + i with max{0, ^ — |} < e < ^ — p G (2, 6) and 
some g e [2, 6], that 



J{t)\u\^u ( , y s^C\\u\\l„ \\J(t)u\\ 



LP 



From the Sobolev embedding theorem and the Holder inequality, it yields 



J(t) lulu 



L^(i-p«) (0,T;l(i-p=) ) 

\\mu\\ 

Similarly, we have 

H{t) |u|^M 



L"'(t^) (0,T;l(t^) ) 



and 



|u|^ u 



For convenience, we denote 

III"IIIg - IIwIIg + II-^WwIIg + II^W^IIg' 

where G denotes a normalized space. Thus, we have 



|m|^ u 



<CT^ 2^ -.1) ||u||?^,,^.„i, iiiMiii 



(4.1) 



l^(t^) (o,T;l(t^) 

LL°'(0,T;ffi) lir"lllL^(p)(0,T;LP) 

For any p G (2, 6) and M ^ 2C Ijuolls' define the workspace (P, d) as 
V := {u: |||u||lioo(Q jn.^2)(-,^",(p)(o,T;Lp) ^ ^^}' 

with the distance 



= Ilk - w|||i^(p)(o_-r.ip) . 

It is clear that (I?, d) is a Banach space. Let us consider the mapping ,^ : {V, d) —> {V, d) 
defined by 

5^ : u{t) ^ S{t)u„ -iP f S{t)S*{s) \u\^ u{s)ds. 



For u e (X*, d), by the commutation relation fTTH)) . ([TTT)) . ([m)) . Proposition 
and the nonlinear estimate (|4.ip . we obtain 



III"^w||IlMp)(0,T;Lp) «^C'll"o|ls + ll"llL~(0,T;ffi) III"! 

(4.2) 

+ \\Oj{t)u\\ + \\OH{t)u\\ 

<M/2 + (CT^"5^^^7V//2 _^ 2CT)M 

<M, (4.3) 
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where we have taken T e (0,7r/4tj] so small that CT y'p^ IVP + 2CT ^ 1/2. Similar 
to the above, a straightforward computation shows that it holds 



1 



L=°(0,T;H1 



lL°°(0,T;i/l) 



2CT 



\\U — V\ 



Liip){0,T;LP) 



CT'-'!'- — M^ + 2CT 



d{u, v) 



^^d{u,v). 



(4.4) 



Hence, ^ is a. contracted mapping from the Banach space {V, d) to itself. By the Banach 
contraction mapping principle, we know that there exists a unique solution u e (I?, d) 
to (|1.5p - p.6p . In view of the conservation laws stated in Proposition l3.H we can use the 
standard argument to extend it uniquely to a solution at the interval [0, 7r/4a;] which 
satisfies, for any t E [0,7r/4ti>] and p e (2,6), that 



u{t, x), J{t)u{t, x), H{t)u{t, x) £ C(0, tt/Alu; L'^{M.^)) n L'^^p^O, tt/Auj; LP{W^)). 



Then, we can extend the above solution to a global one by translation. In fact, 
in order to get the solution in the interval (7r/4cL;, 7r/2w], we can apply a translation 
transformation with respect to the time variable t such that the initial data u^tt/Auj) are 
replaced by u{0). Let u{t,x) := u{t— 7r/4a;,a;), then we have from the original equation 
with initial data u(tt/Auj,x) 

iut + -Ail = — \x\'^ U + /3 u - ujL^u, x £ M^ t ^ 0, (4.5) 
u{0, x) = uo{x) := ii(7r/4tj, a;), x G M^. (4.6) 

In the same way, we can get a solution u{t,x) of ()4.5P " (|4.6p for t G [0,7r/4w]. It is 
also a solution u{t,x) to ()1.5|) " (ll.6|) for t G [n / Auj , tt / 2uj] and it is unique. Thus, by an 
induction argument with the help of those conserved identities stated in Proposition l3.H 
we can obtain a global solution u{t, x) to (|1.5p - (|1.6p satisfying for any T G (0, oo) 

u{t, x), J{t)u(t, x),H{t)u{t, x) e C(M; L^{W^)) n L^'-p^ (0, T; LP{R^)). 

Therefore, we have completed the proof of the main theorem. 
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